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Further Evidence for an Elliptical Instability in Rotating Fluid 

Bars and Ellipsoidal Stars 

Shangli Ou^'^ Joel E. Tohline^ & Patrick M. Motl^ 
ABSTRACT 

Using a three-dimensional nonlinear hydrodynamic code, we examine the dy- 
namical stability of more than twenty self-gravitating, compressible, ellipsoidal 
fluid configurations that initially have the same velocity structure as Riemann 
S-type ellipsoids. Our focus is on "adjoint" configurations, in which internal fluid 
motions dominate over the collective spin of the ellipsoidal figure; Dedekind-like 
configurations are among this group. We find that, although some models are 
stable and some are moderately unstable, the majority are violently unstable 
toward the development of m = 1, m = 3, and higher-order azimuthal distor- 
tions that destroy the coherent, m = 2 bar-like structure of the initial ellipsoidal 
configuration on a dynamical time scale. 

The parameter regime over which our models are found to be unstable gen- 
erally corresponds with the regime over which incompressible Riemann S-type 
ellipsoids have been found to be susceptible to an elliptical strain instability 



(ILebovitz fc Lifschitzlll996l ). We therefore suspect that an elliptical instability is 
responsible for the destruction of our compressible analogs of Riemann ellipsoids. 
The existence of the elliptical instability raises concerns regarding the final fate 
of neutron stars that encounter the secular bar-mode instability and regarding 
the spectrum of gravitational waves that will be radiated from such systems. 

Subject headings: Riemann S-type ellipsoids — Galaxies: Rotating bars — Star 
Formation: Rotating ellipsoids — Neutron stars — Elliptical fluid instability — 
computational astrophysics 
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Introduction 



In a previous paper (IOull2006l . hereafter paper I), we introduced a new self-consistent- 
field technique that is capable of constructing three-dimensional (3D) models of incom- 
pressible Riemann S-type ellipsoids and compressible triaxial configurations that share the 
same velocity fields as those of Riemann S-type ellipsoids. These compressible triaxial con- 
figurations represent fairly good quasi-equilibrium states and can be used to examine the 
dynamical stability of Riemann S-type ellipsoids in the nonlinear regime. With this in 
mind, a subset of our compre ssible models have been evolved in a 3D hydrodynamics code 
( iMotl. Tohline. fc Frankll2002l ) for 20 — 40 dynamical times to test their dynamical stability. 
In this paper, we present results from these hydrodynamic simulations. 

Before going into the details of our simulations, it is worthwhile spending some time 
reviewing previous studies of the stability properties of Riemann S-type ellipsoids. Chan- 
drasekhar and Lebovitz carried out a vigorous analysis of the stability of Riemann S-type 
ell ipsoids with respect to second and third order harmonics (these results are summarized 



m 



Chandrasekhan (119691 )). Some configurations, including both rigidly- rotating Jacobi and 



stationary Dedekind sequences, were found to be subject to a pear-mode (azimuthal mode 
number m = 3) instability beyond a certain point while proceeding away from the axisym- 
metric Maclaurin spheroid sequence. Part of the Jacobi and Dedekind sequences were also 
found to be unstable to a dumbbell-shaped (azimuthal mode number m = 4) mode; this has 
provided some theoretica l foundation for the fission the ory of the formation of binary stars 
jDurisen fc Tohlinel Egssi : llebovitzl Il987l : IXohhnel I2OO2I ) . These studies have assumed that 
the evolution of a rotating gas cloud toward the bifurcation points of the m = 3 pear-mode 
or m = 4 dumbbell-mode proceeds in a quasi- stationary manner and on a time scale that is 
fairly long compared to the dynamical time scale. 



However, more recent linear stability studies (ILebovitz fc Lifschitzlll996l : iLebovitz fc Saldanha 



19991 ) have revealed a new elliptical strain instability in Riemann S-type ellipsoids that is 
associated with the noncircular fluid streamlines that describe the internal motion of these 
configurations. According to these studies, a fairly large fraction of all rotating ellipsoids 
appear to be susceptible to the development of this elliptical instabihty. It has a growth 
rate that is slow compared to the rotation rate but is nevertheless dynamical. The fastest 
growth rates are found among the adjoint configurations, in which the internal motion dom- 
inates over the collective patt ern motion of the rot ating bar. This el liptical instability was 
first discovered and studied by iPierrehumben ( 



dimen sional (2D), incompressible flu id flows. 



19861) andlBavlyi (119861) in the context of two- 
Goodman (jl993 ). Lubow. Pringle fc Kerswell 



(I1993I ) , and iRyu fc GoodmanI (Il994j ) extended it to astrophysics to study tidally distorted 
accretion disks. They found that the instability is three-dimensional (depends on the vertical 
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dimension) and approximately incompressible, and that it might also give birth to m = 1 
inter nal waves in the absence of viscosity. It is subsequent to this work on accretion disks 
that iLebovitz fc Lifschitzl (119961 ) conducted their study of Riemann S-type ellipsoids, for 
which elliptical flows are generic. 

The discovery of an elliptical strain instabil ity in rotating ellipsoidal (bar-like) structures 
brings into question ( Lebovitz fc Lifschit j|l996] ) whether the paradigm for evolution of a sec- 
ularly unstable star drive n by gravitational-radiation reaction (GRR) forces is viable. Ac- 
cordi ng to this paradigm (I Chandr asekharl 1 1 9 70 : iDetweiler fc LindblomI 119771 : iLai fc Shapiro 
19951 ). such a star - in the form of a Riemann S-type ellipsoid - would evolve on a secular 
time scale toward the Dedekind sequence (for which the pattern frequency of the bar-mode 
a; = 0), at which point the evolution would stop because the gravitational field no longer 
exhibits a time- varying quadrupole moment. The analysis of iLebovitz fc Lifschita (119961 ) 
shows that the parameter space occupied by the elliptical instability encompasses most of 
the Dedekind sequence. Because the elliptical instability develops on a dynamical time scale, 
it might prevent a GRR-driven figure from proceeding further to the Dedekind sequence, es- 
pecially if th e elliptical instability develo ps into the nonlinear regime. Because the results 
presented by ILebovitz fc Lifschitzl (119961 ) are based on linear analysis, they were not sure 
what the nonlinear outcome of the elliptical instability would be. 

In a recent nonlinear study of the development of the secular bar-mode instability by 
Ou. Tohline. fc LindblomI (120041 ). an initially uniformly rotating neutron star was driven by 
GRR forces into a differentially rotating bar-like configuration that had a very low pat- 
tern frequency. This bar-like configuration appeared to be a compressible analogue of an 
adjoint configuration among the Riemann S-type ellipsoids. However, the coherent bar- 
like structure was destroyed very quickly by the nonlinear development of an unexpected 
high- order instability that generated "turbulence" throughout the internal structure of the 
bar. lOu. Tohline. fc LindblomI (120041 ) suggested that this turbulent-like instabihty was re- 
lated to the elliptical instability, but no definite conclusion could be drawn because only 
one model was studied. It was impractical to pursue a full investigation across the entire 
parameter space of possible ellipsoidal models because, at the time, the generation of one 
single equilibrium bar-like model was computationally expensive; the construction of even 
one equilibrium ellipsoidal model required following the GRR-driven evolution of an initially 
axisyr nmetric model, which takes a v e ry long time even with an artificially enhanced GRR 
force JOu. Tohline. fc LindblomI booi Ishibata fc Karinoll2004h . 



With the new self-consistent-field technique introduced in Paper I, we have been able 
to quickly build a large number of nearly incompressible, quasi-equilibrium triaxial models. 
These compressible models contain the same internal velocity fiow-fields as incompressible 
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Riemann ellipsoids. And although they are not exact steady-state configurations, these mod- 
els provide reasonably good initial states for a study of the elliptical instability because the 
elliptical instability develops on a dynamical time scale. In what follows, we present hy- 
drodynamic simulations of a variety of these nearly incompressible, rotating bar-like models 
to test their dynamical stability. We find that a significant fraction are unstable to what 
appears to be the elliptical instability. We introduce our initial models and analysis method 
in §2; we present results from hydrodynamic evolutions for different configurations in §3; and 
in §4, we conclude and discuss issues relevant to the evolutionary path of secularly unstable 
neutron stars. 



2. Initial Models and Methods 

2.1. Initial Models 

We follow the notation defined in Paper 1, where a, b, and c are the three principal 
semiaxes of an ellipsoid and c is the rotation axis. In a frame of reference whose origin 
coincides with the center of the ellipsoid and that is rotating with the pattern frequency of 
the ellipsoid, u, the velocity field of a Riemann S-type ellipsoid takes the form, 

V = X{ay/b,—bx/a,0) , (1) 

where A is a constant that determines the degree of internal motion of the fiuid. This 
velocity field is designed so that velocity vectors are everywhere tangent to a set of self- 
similar concentric ellipses. (In the uniform-density, incompressible Riemann configurations, 
these concentric ellipses also identify equipotential contours.) When \uj\ > |A|, a model is 
said to be a "direct" model, in which the figure rotation dominates over internal motions. 
Conversely, when \uj\ < |A|, the configuration represents an "adjoint" model, in which the 
internal fiuid motion dominates over the figure rotation. 

In Paper I, we constructed a variety of compressible models that have the same velocity 
field as is found in incompressible Riemann S-type ellipsoids, that is, with the velocity field 
defined by Eq. ([1]). These models satisfied the steady-state Euler's equation exactly, but 
only satisfied the steady-state continuity equation approximately. However, as discussed in 
Paper I, it was found that the deviation from the steady-state continuity equation was small 
for models that obeyed a relatively stiff equation of state (EOS). For our present study, we 
have used the same self-consistent-field technique described in Paper I to construct more 
than forty quasi-equilibrium, compressible analogues of Riemann S-type ellipsoids to serve 
as initial states for a series of hydrodjTiamic simulations. The models were constructed 
using two different EOS {n = 0.5 and n = 1 polytropes), with a variety of initial axis 
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ratios - Figure [H] displays the range of selected (6/a, c/a) pairs - and, for a given set of 
axes, both "direct" and "adjoint" configurations were constructed. This set of models was 
sel ected to cover a reas o nably large fraction of the parameter space that was investigated 
by iLebovitz &: Lifschitzl (119961 ) , keeping in mind that their linear analysis was confined to 
incompressible figures. 

The stability of roughly twenty-five of these initial models was examined using a 3D 
hydrodynamic code, as explained below. We focused on analyzing the stability of adjoint 
models with an n = 0.5 polytropic EOS, but several direct models with an n = 0.5 EOS and 
one adjoint model with an n = 1 EOS were included in the mix. For the sake of brevity, 
only the results of six model evolutions are presented in detail in ^ Table [1] lists the 
key parameters for these six representative models. In addition to the name that has been 
assigned to each model and the values of five parameters that already have been defined 
(n, 6/a, c/a, cj. A), Table [1] specifies each initial model's ratio of rotational to gravitational 
potential energy T/|l^|, total angular momentum Jtot, mean density pmean , and the ratio of 
the pattern period Tp = 27r/co' to the dynamical time Tdyn = [TrGpmcan]^^^^- The measured 
growth rates Ti and T3 of unstable modes with azimuthal mode numbers m = 1 and m = 3 
are also tabulated (see §3 for further elaboration). Dimensional quantities are in the hydro 
code unit for which the gravitational constant G, the central density pc-, and the cylindrical 
radius of the entire grids Wgrid are all set to unity. If the first letter in a model's name is 
D, the model is a "direct" configuration, whereas an A denotes an "adjoint" configuration. 
Among this group of six models, we note that only model A067 has a polytropic index 
n = 1\ the other models all have n = 0.5. 



Table 1. Parameters of selected initial models 



model 


n 


h/a 


c/a 


UJ 


A 


T/\W\ 


^tot 


Pmcan 




^l/^dyn 




D105 


0.5 


0.59 


0.487 


0.911 


0.114 


0.105 


0.00809 


0.526 


1.4 






D080 


0.5 


0.49 


0.487 


0.885 


0.230 


0.080 


0.00497 


0.529 


1.5 






AOlO 


0.5 


0.74 


0.821 


-0.610 


-0.875 


0.010 


7.24e-3 


0.543 


2.1 






A067 


1.0 


0.90 


0.692 


-0.168 


-0.796 


0.067 


0.0087 


0.280 


5.6 






A134 


0.5 


0.74 


0.487 


-8.1e-3 


-0.911 


0.134 


0.0133 


0.519 


158 


3.3 


3.4 


AlOO 


0.5 


0.59 


0.487 


-0.114 


-0.911 


0.100 


7.21e-3 


0.526 


11 


1.4 


1.3 
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2.2. Methods 

In order to examine their relative stabihty, each model was evolved on a cylindrical coor- 
dinate mesh using the 3D, Newtonian, finite- vo l ume c omputational fiuid dynamics technique 
described in detail by lMotl. Tohline. fc FrankI (120021 ). In each simulation, the adopted grid 
resolution was 66 x 102 x 128 in the w (cylindrical radius), z (vertical), and (azimuthal) 
directions, respectively. 

We monitor the time-evolution of the general structure of each model by measuring the 
Fourier amplitude of various "modes" having azimuthal quantum numbers m in the following 
fashion: At each instance of time, the azimuthal density distribution in a ring of fixed R and 
z can be decomposed into a series of azimuthal Fourier components via the relation. 



p(i?,^,0)= Cm{R.z)e' 



iTn<f> 



(2) 



m=—oo 



where the complex amplitudes Cm are defined by the expression. 



1 



2n 



Cm{R.z) = — / p{R,z,<i))e-'^^d(t). 
2vr Jo 



(3) 



In our simulations, the time-dependent behavior of the magnitude of this coefficient, \Cm\i is 
monitored at a variety of (i?, z) locations to measure the growth-rate of various modes. The 
time- varying amplitudes shown in our figures, below, are averages of \Cm\ over the entire 
volume. 

Although this diagnostic tool permits us to study the behavior of structure having a 
wide variety of m values, in practice we have focused our attention on the lowest order 
modes for this study. For example, a trace of the time-evolution of IC2I lets us monitor the 
overall ellipticity of each model. The m = 2 mode should dominate initially because of each 
model's initial bar-like structure. If this amplitude remains at a fairly constant level in a 
given simulation, we conclude that the configuration is stable; otherwise, it is unstable. In 
practice, we have found that when the m = 2 amplitude drops significantly in a model, that 
model simultaneously displays growth of other azimuthal modes. 



In their analysis of the elliptical instability, iLebovitz fc Lifschitzl (119961 ) generally found 
that the unstable mode with the fastest growth rate was the m = 3 mode. In their study, 
the m = 1 mode was not discussed because, in incompressible configurations, an m = 1 
mode normally refers to a translation of the center of mass of the configura tion. However, in 
the study of an elliptical instabi li ty in tidally distorted accretion disks by iGoodmanI (119931 ) 
and iLubow. Pringle fc Kerswelll (119931 ). it was found that an eccentric m = 1 mode may 
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be important. Further investigation shows that the instabihty has sor ne dependence on 
the t hird (vertical) dimension, because it disappeared in 2D simulations (IRyu fc Goodman 



1994J ). With these results in mind, we have focused our attention on the development of the 



m = 1 and 3 modes in our 3D ellipsoidal model evolutions. 



3. Results From Hydrodynamic Simulations 



We present our results from hydrodynamic evolutions for direct and adjoint configura- 
tions in the next two subsections, respectively. 



3.1. Results For Direct Configurations 



In this subsection, we present results from two direct configurations, models D105 and 
D080, which were evolved for about 24 Tdyn and 17 Tdyn, respectively. Figure [1] shows 
snapshots in time of equatorial-plane iso-density contours and the inertial-frame velocity 
field for model D105. (In the online version of this paper, the figure caption includes a link 
to an MPEG movie showing the entire evolution of this model.) This model is stable; the 
coherent bar-like structure and elliptical flows are well maintained throughout the 24 Tdyn 
followed by our simulation. The time-varying Fourier mode amplitudes displayed in Figure 
|2] show that the m = 2 bar-mode remains the dominant mode throughout this evolution, 
with an amplitude that remains approximately three orders of magnitude higher than those 
of the m = 1 and 3 modes. Although the m = 4 mode has a larger amplitude than the 
odd modes, the fact that its behavior follows that of the m = 2 mode suggests that it is a 
harmonic of the m = 2 mode and that there is no independent m = 4 mode. 

Although we have categorized model D105 as stable, a slow decay of the bar-mode 
amplitude is evident in Figure 2. Because other modes, including the m = 1 and 3 modes, 
do not grow throughout the evolution, this slow decay of the bar-mode probably results from 
the small violation of the steady-state continuity equation in our initial model configuration. 
(This, in turn, results from the fact that the 3D iso-density surfaces in our compressible 
models deviate from self-similar ellipses or homeoidal shells; see Paper I for a more detailed 
discussion.) Small violations in the surface layers cause very low density material to be shed 
outward from both ends of the object and form two low density spiral waves (see the panel 
labeled t = 5 in Figure 1). Because this material carries relatively high specific angular 
momentum, the overall system loses angular momentum, but only on a very long time scale. 
We suspect that this is the major cause of the gradual decay of the m = 2 bar-mode in this 
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model evolution. 

Figure [3] shows a similar Fourier mode analysis for model D080. Again, the m = 2 
bar-mode dominates the entire evolution and shows a slow decay; amplitudes of the odd 
azimuthal modes remain at roughly constant "noise" levels. 

The results from the evolution of these two direct configurations suggest that, although 
our initial model configurations do not satisfy the steady-state continuity equation exactly, 
the models are still fairly good quasi-equilibrium states. In both cases, the coherent bar-like 
structures were well maintained through ~ 20 Tdyn- Hence, we are confident that initial 
models generated by the technique described in Paper I can be used to study the dynamical 
stability of compressible Riemann S-type ellipsoids. Furthermore, the fact that the elliptical 
instability sets in on a dynamical time scale will allow us to use these models to detect the 
elliptical instability, if it exists. 



According to iLebovitz fc Lifschita (119961 . see their Figure 2 for a summary of their 



analysis of direct configurations), the elliptical instability only arises in a small portion 
of the parameter space that is occupied by direct configurations, and its growth rate is 
relatively slow there. This makes it relatively difficult for our current study to identify 
the exact parameter-space boundaries between stable and unstable direct configurations. 
Furthermore, the method described in Paper I can only be used to construct compressible 
models having axis ratios b/a > 0.5, whereas most of the domain that we expect to be 
effected by the elliptical instability in direct configurations has b/a < 0.5. This also makes 
direct configurations unsuitable for identifying the existence of the elliptical instability. 

Therefore, we have chosen to focus our study on an analy sis of adjoint configurations 



where the domain of the instability is much wider, according to ILebovitz fc Lifschitzl (Il996l . 
see especially their Figure 3a, which is reproduced here as Figure [T3i) . More importantly, 
the growth rate of the instability is expected to be faster in adjoint configurations. This 
helps us to conduct such a survey with a limited evolution time that is allowed by our 
current computational resources. In the next subsection, we present simulations of adjoint 
configurations, which indeed reveal the development of turbulence-like fiows that we suspect 
are associated with the elliptical instability. 



3.2. Results of Adjoint Configurations and Existence of a Turbulent-like 

Instability 



In this subsection, we present results from 3D hydrodynar nic evolutions of the four ad- 

joint configurations listed in Tabledl In light of previous studies jGoodmaJl993l :l ILebovitz Sz Lifschitzl 
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19961 ). it is expected that the m = 1 and 3 modes will be the fastest growing modes, if an el- 
liptical instability sets in. Therefore, rapid growth of these two modes can serve as a compass 
pointing to the elliptical instability. With this in mind, in the analysis below we will catego- 
rize our results based on the behavior of these two modes. Based on the final outcome and 
behavior of the m = 1 and 3 modes, evolutions of our initially quasi-equilibrium compress- 
ible adjoint configurations can be divided into three categories: stable models, moderately 
unstable models, and violently unstable models. 

The first category — stable models — consists of evolutions in which the internal el- 
liptical flow remains stable. One evolution in this family is model AOlO, with b/a = 0.74. 
Note from the information in Table [1] that this model has a somewhat prolate structure 
with c > b. Figure H] shows 2D, equatorial-plane iso-density contours of this model at four 
different times. (The figure caption in the online version includes a link to an MPEG movie 
showing the entire evolution of this model.) The ellipsoidal shape is well preserved, at least 
up through the end of our simulation ~ 26 T^yn, and the internal elliptical flow pattern 
remains stable. As suggested by the Fourier mode analysis shown in Figure [5l the m = 2 
bar-mode dominates throughout the entire evolution. There is no indication that the m = 1 
and 3 modes are growing; their amplitudes remain fairly constant at "noise" level throughout 
the evolution. Figures [6] and [7] show similar plots for another stable model, namely, model 
A067 with b/a = 0.90. This model was perturbed with an initial low-amplitude m = 3 
perturbation, but its bar-like structure remains stable during the evolution. We note that 
the bar-mode appears to undergo a low-amplitude oscillation. As we have discussed in Paper 
I, initially small deviations from the steady-state continuity equation can generate surface 
waves in adjoint models. The oscillations observed here in the amplitude of the m = 2 bar- 
mode is indeed a reflection of these low-amplitude surface waves traveling along the elliptical 
surface. The initial violation of the steady-state continuity equation is larger in models with 
a softer EOS (see Paper I). It is therefore not surprising that oscillations of the bar- mode 
amplitude in model A067 are somewhat larger than in other models because model A067 
was constructed using a relatively soft {n = 1 polytropic) EOS. 

The second category of adjoint evolutions consists of configurations that are moderately 
unstable. In these models, the elliptical flow appears to remain stable as judged by the 
behavior of the m = 2 Fourier mode. However, the m = 1 and 3 modes are observed to 
grow on a long time scale. Model A134 with b/a = 0.74 belongs to this family. Figure [8] 
shows 2D iso-density contours of this model in the equatorial plane at different evolutionary 
times. (An MPEG movie showing part of the evolution of this model is available in the figure 
caption of the online version.) We note that this model has an almost vanishingly small bar 
pattern frequency, uj = —0.0081, so it is very nearly a Dedekind-like object; as seen in Figure 
[HI the bar-like pattern shifts retrograde (clockwise) through only ^ 7r/4 radians over 15 Tdyn! 
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Through this initial ~ 15 Tdyn, the coherent elhpsoidal structure and eUiptical flow pattern 
remain stable. However, as revealed by the Fourier mode analysis in Figure [9l at later times 
the m = 1 and 3 modes both exhibit exponential growth, albeit at a very slow growth rate. 
We note that the amplitudes of these odd modes appear to be tied to each other; even their 
individual growth rates are almost the same. (As recorded in Table [H the e-folding time for 
these odd modes are Ti = 3.3Tdyn, and T3 = 3.4Tdyn.) This suggests that the mechanism 
that is responsible for driving these modes is probably the same. Their amplitudes do not 
exceed that of the m = 2 mode at the end of this simulation. Hence, we have identified this 
model as being moderately unstable. Across the parameter space of our investigation (see, 
for example. Figure [T^ . models belonging to this category can be regarded as a transition 
from stable to violently unstable models. 

The third category consists of models that are violently unstable to an instability that 
leads to apparently turbulent flow throughout the configuration. The m = 2 bar-mode is 
quickly destroyed in a few dynamical times; simultaneously, both the m = 1 and 3 modes 
exhibit extremely fast exponential growth. One of the most unstable models in our investi- 
gation is AlOO with b/a = 0.59. The time-evolution of equatorial-plane iso-density contours 
of this model is shown in Figure [TDl (The entire evolution is shown in a MPEG movie in 
the figure caption of this plot in the online version of this paper.) The initial ellipsoidal 
structure exists for only ~ 3Tdyn, then it is suddenly destroyed by some type of "turbulent" 
instability. After the bar-like structure is destroyed, the whole configuration settles into an 
entirely turbulent state, exhibiting small-scale vortices and eddies. The Fourier mode anal- 
ysis shown in Figure [TT] looks distinctly different from all previous mode plots: the m = 2 
mode decays very quickly in parallel with rapidly growing m = 1 and 3 modes, which finally 
surpass the m = 2 mode in amplitude. These distinct features suggest that the initially el- 
lipsoidal structure of model AlOO was destroyed by a process that is entirely different from 
the gradual decay of the bar-mode amplitude seen in our models that have been categorized 
as stable. As in model A134, the growth time for both the m = 1 and 3 modes are very close 
to each other (as recorded in Tabled! Ti = 1.4Tdyn and T3 = 1.3Tdyn), but they ultimately 
saturated at different levels. We evolved model AlOO for a long time (~ 50Tdyn) until the 
system reached a new, nearly axisjTumetric state. 

The turbulent behavior of model AlOO is very similar to the outcome of an instabil- 
ity tha t arose during the late evolutio n of the GRR-driven neutron star from a previous 



study (jOu. Tohline. fc Lindblomll2004l ). In an effort to compare this earlier simulation re- 
sult with our present one, we have plotted in Figure [T2] a power spectrum of the various 
azimuthal modes for model AlOO at time t = and t = 45Tdyn in the evolution. This 
power spectrum looks su r prisiii gly similar to the power spectrum that was generated by 



Qu. Tohline. fc Lindblonj (120041 ) late in their model's evolution. The spectrum is dominated 
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by even modes in the beginning as a result of the elhpsoidal (bar-hke) structure; whereas 
a turbulence-hke power law distribution dominates at late times, after the energy initially 
stored in the m = 2 bar-mode has cascaded into all other Fourier modes. This strongly sug- 
gest s that the violent instability seen in our present study is the same instability that arose in 
the lOu. Tohline. fc LindblomI (120041 ) investig ation. Because this turbu l ence-l ike instability 



has arisen in two quite different contexts — lOu. Tohline. fc LindblomI (120041 ) discovered it 



while studying the non-linear evolution of an initially axisymmetric star that was driven to 
a bar-like configuration by GRR forces, whereas we have discovered it during purely hydro- 
dynamic evolutions of compressible triaxial models that have been purposely constructed 
to resemble adjoint configurations of Riemann S-type ellipsoids — it seems likely that the 
instability is generic to adjoint configurations. 

In an effort to understand whether or not this insta bility occupies the s a me re gions 
of parameter space as the elliptical instability studied by iLebovitz fc Lifschita (Il996l ). our 
hydrodynamic investigation has included models that cover as much parameter space as 
possible. In Figure [TU we have mapped the results from 20 of our adjoint model evolu- 
tions onto the (6/a, c/a) plane. (For the sake of brevity, a detailed discussion of only 6 of 
these model evolutions has been presented, above.) Models found to be violently unstable 
are marked by asterisks, stable models are marked by diamonds, and moderately unstable 
models are marked b y plu s signs. This figure should be compared with Figure 3a from 
Lebovitz fc Lifschitzl ( 119961 ) - reproduced here as Figure [13] - which shows the instability 
domain occupied by the m = 3 mode in the adjoint family of incompressible, Riemann S- 
type ellipsoids. It is very e ncouraging that our regions of instability are consistent with the 
Lebovitz fc Lifschitzl (119961 ) diagram. More specifically, both investigations find that unsta- 
ble configurations occupy a large portion of parameter space above the Dedekind sequence 
when b/a < 0.9. The growth time also appears to be related to the axis ratio, b/a: mod- 
els with smaller b/a undergo a more violent instability. Despite the fact that the method 
described in Paper I prevents us from constructing compressible models with c/a < 0.4 or 
b/a < 0.5, our results already indicate that models that are most susceptible to the instabil- 
ity occupy a smaller fraction of the available parameter space below the Dedekind sequence 
than above it. The two models that we have categorized as being moderately unstable (the 
plus signs in Figure [T^ mark the transition from stable to violently unstable domains. In 
particular, the fact that model A134 with {b/a, c/a) = (0.74,0.487) is only moderately 
unstable suggests that we may be seeing evidence for the "tongue" structure displayed in 
Figure [13] that covers most of the Dedekind sequence. As a counter e xample, however, we 
note t hat all of our models with b/a = 0.89 appear to be stable whereas [Lebovitz fc Lifschitz 
( 1l996[ ) found that some models with b/a 0.9 fall in the instability domain covered by a set 
of small extended tongues. This apparent discrepancy probably arises because the growth 
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rate of the instability is sufficiently slow in this regime that we are unable to identify the 
instability, given the limited evolutionary time of our simulations. 

Because our models are constructed on a rather coarse, discrete parameter grid, we are 
unable to map out the instability domain with the type of fine structure shown in Figure [T3l 
However, the fact that the instability domain identified in our present s tudy and displayed in 



Figure [Ml by and large agrees with the instability domain identified by iLebovitz &: Lifschitz 



(119961 ) gives us confidence that the turbulence- like behavior observed in our violently unstable 



models arises from the elliptical instability. 



4. Conclusions and Discussions 



We have conducted a survey of the nonlinear hydrodynamic stability of compressible 
triaxial configurations that are initially in quasi-equilibrium and share the same velocity 
field as that of Riemann S-type ellipsoids. Our simulations show that many of these models 
are indeed very good long-hved, quasi-equilibrium states. In addition, we have found that 
a turbulence-like dynamical instability arises over a fairly large region of the examined pa- 
rameter domain (see Figure IT^ . This instability is probably the same instability observed 
i n the late evolution of a neutron star that was driven to a bar-like structure by GRR forces 
(jOu. Tohline. &: Lindblomll2004l ). The characteristics exhibited by this instability and the 
domain it occupies in the examined parameter space agree qualitatively with the elliptical 
instab ility that has been discove red in previous linear studies of incompressible configura- 
tions (iLebovitz fc Lifschitzlll996l ). Therefore, we suspect that this instability is the elliptical 
instability that appears to develop generically in fluid flows with elliptical stream lines. 

As a result of this instability, we found that all odd azimuthal modes grew on a dynamical 
time scale, but were dominated by the m = 1 and 3 modes, whose amplitudes eventually 
surpassed the m = 2 bar-mode and destroyed the ellipsoidal structure of each violently 
unstable model. It appears that the growth rate for unstable modes is higher in models with 
smaller b/a values. This makes sense because the elliptical instability ought to disappear 
in circular flows with b/a = 1. In a li n ear an alysis of the elliptical instability in Riemann 
S-type ellipsoids, iLebovitz fc Lifschita (ll996l ) did not discuss the behavior of t he m = 1 
mode. However, the appeara nce of an eccentric m = 1 mode in accretion disks ( GoodmanI 
19931 : iRyu fc Goodmanlll994l ) is consistent with our discovery here that the m = 1 mode can 
also develop in compressible flows and drain energy from the m = 2 mode. 



The existence of the elliptical instability raises concerns regarding the final fate of neu- 
tron stars that encounter the GRR-induced secular bar-mode instability. According to pre- 
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vious theoretical investigations (jPetweiler fc Lindbloml 119771 : iLai fc Shapiro! Il995l ). such a 
star would be expected to evolve through a sequence of Riemann S-type ellipsoids toward 
the Dedekind sequence that has a vanishingly small pattern frequency as viewed from the 
inert ial reference frame. However, iLebovitz &: Lifschitz Jl996h have shown that most of the 
Dedekind sequence falls in the parameter regime where fluid configurations should be sus- 
ceptible to the elliptical instability; if any star enters this region of parameter space, the 
elliptical instability would se t in on a dynamical time scale. It was unclear in the linear 
study of iLebovitz fc Lifschitzl (119961 ) whether the elliptical instability would be mild or vio- 
lent in the nonlinear regime. This picture was partially answered by the nonlinear evolution 
of one rota ting neutron star that was sus ceptible to the GRR-induced secular bar-mode 
instability (lOu. Tohline. &: Lindbloml 120041 ) : the bar-like configuration was observed to be 
destroyed very quickly when the pattern frequency of the bar dropped to nearly zero. The 
evolutions of compressible triaxial models presented in this paper further clarify what the 
nonlinear outcome will be of the elliptical instability. We have found that the instability can 
be violent and capable of destroying the ellipsoidal structure of the star. Therefore, it seems 
unlikely that a secularly unstable star driven by GRR forces will actually be able to evolve 
along a sequence of Riemann S-type ellipsoids toward a Dedekind-like configuration. 



At the final stage of the evolution reported by lOu. Tohline. &: Lindbloml (120041 ). after 
the bar-like structure was destroyed and the star returned back to a nearly axisymmetric 
state with strong differential rotation, its r/|W^| value was still above the critical limit 
of 0.14 for the secular bar- mode instability. Hence, it is possible that such a star would 
again attempt to evolve toward the Dedekind-sequence under the influence of GRR forces. 
Cycles of instability could continue until the star's rotational energy is drained below the 
critical limit. If this picture is correct, it would provide a very interesting, recurrent source of 
gravitational-wave radiation for ground-based gravitational- wave detectors such as the Laser 
Interferometer Gravitational- wave Observatory (LIGO) and its sister instruments worldwide. 
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Fig. 1. — Equatorial- plane iso-density contours and velocity fields in the inertial frame for 
model D105. Density contours correspond to p/pmax = 0.01, 0.1, 0.5, and 0.8 from the out- 
ermost to the innermost shell, and time is given in units of T^yn- The pattern rotation of the 
overall ellipsoidal configuration and the internal fiuid motion are both prograde (counter- 
clockwise) in this "direct" model configuration. An accompanying MPEG movie shows the 
entire evolution of this model, through t ^ 24 T^yn; 3D, rather than 2D, iso-density contours 
are shown in the movie without the accompanying fiow field. 
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Fig. 2. — Time-evolution of the amplitudes \Cm\ of the m — 1,2, 3, and 4 modes for model 
D105. The amphtudes of the m — 2 and 4 modes decay gradually; the m — 1 and 3 modes 
remain at a fairly constant, low-amplitude level, which suggests that this "direct" model is 
stable against the elliptical instabihty. 
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Fig. 3. — The same as Figure [21 but for model D080. This "direct" model also appears to 
be stable against the elliptical instability. 
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Fig. 4. — The same as Figured], but for model AOIO. In this "adjoint" model, the internal 
fluid motion is prograde (counterclockwise), whereas the overall motion of the ellipsoidal 
configuration is retrograde. An accompanying MPEG movie shows the entire evolution of 
this model, through t ^ 26 T^yn- 
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Fig. 5. — The same as Figure [2], but for model AOIO. All the modes appear to be stable in 
this "adjoint" model. 
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Fig. 6. — The same as Figure [H but for model A067. In this "adjoint" configuration, the 
fiuid is moving prograde (counterclockwise) while the overall ellipsoidal pattern is spinning 
retrograde. 




Fig. 7. — The same as Figure [21 but for model A067. All the modes appear to be stable. 
Note that a low- amplitude, m = 3 perturbation was imposed on this model at time t = 0. 
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Fig. 8. — The same as Figured! but for model A134. The internal fluid motion is prograde 
(counterclockwise), but this "Dedekind-like" model has an almost vanishing angular pattern 
frequency, uj = —0.0081; as seen here, the figure pattern shifts clockwise by only ^ 7r/4 
radians in the first 15 T^yn- From the four snapshots displayed here, the bar structure seems 
to be very stable, but from the mode analysis shown in Figure [HI we conclude that the model 
is moderately unstable. An accompanying MPEG movie shows the entire evolution of this 
model, through t ^ 40 T^yn- 
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Fig. 9. — The same as Figure [2], but for model A134. The amphtude of the m = 2 mode 
remains fairly constant throughout most of this evolution. The m = 1 and 3 modes also 
seem to be stable at early times, but ultimately they both exhibit a period of exponential 
growth (with a relatively slow growth time), which suggests that the model is moderately 
unstable. 
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Fig. 10. — The same as Figure [H but for model AlOO. These frames have been selected 
between 3 < t/T^jn < 8, during which the initial bar-like structure is entirely destroyed by 
some turbulent instability. An accompanying MPEG movie shows the entire evolution of 
this model, through t ^ 50 Tjyn. 
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Fig. 11. — The same as Figure [2], but for model AlOO. The m = 1 and 3 modes grow very 
quickly, in concert with the overall development of turbulent motions inside the configuration. 
At the same time, the m = 2 mode also decays very quickly. The initial bar-like structure is 
destroyed in less than 8 Tdyn- 
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Fig. 12. — The azimuthal Fourier mode power spectrum for model AlOO initially (filled 
circles) and at late times (open circles) . To guide the eye, amplitudes determined for various 
modes at the same time are connected by straight line segments. Initially, all the even modes 
have higher amplitude than their neighboring odd modes, whereas a cascading power-law 
distribution appears at late times. 
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Fig. 13. — Results from Figure 3a of iLebovitz fc Lifschita (119961 ) are redrawn here to pro- 
vide background to our present investigation. Configurations that belong to the adjoint 
family of incompressible, Riemann S-type ellipsoids exist in this (6/a,c/a) parameter do- 
main everywhere between the upper solid curve - drawn from the origin (6/a,c/a) = (0,0) 
to (6/a, c/a) = (1,1) - and the lower solid curve - drawn from the origin to (6/a,c/a) = 
(1,0.303). Dedekind configurations lie along the solid curve that connects the origin to 
(6/a, c/a) = (1,0.583). The jagged, dotted line that connects the upper curve to the lower 
curve marks the boundary between adjoint configurations that are stable (ellipsoids to the 
right of the line) and unstable (ellipsoids to the left of the line) toward the growth of m = 3 



perturbations, as determined from the linear stability analysis of lLebovitz fc Lifschita (119961 ). 
Note that in certain regions, a thin "tongue" of unstable configurations penetrates into the 
region that otherwise contains stable ellipsoids. 
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Fig. 14. — Results from our present investigation are summarized across the {b/a,c/a) 
parameter domain from twenty separate dynamical simulations of compressible, adjoint con- 
figurations. The three solid curves are the same as in Figure [131 Asterisks denote violently 
unstable models, plus signs denote moderately unstable models, and diamonds denote stable 
models. The four models whose evolutions have been discussed in detail in the text can be 
identified by the {b/a,c/a) parameters defining their initial shapes, as documented in Table 
1: AOlO (0.74,0.821); A067 (0.90,0.692); A134 (0.74,0.487); and AlOO (0.59,0.487). Of the 
models identified here, only model A067 was constructed using the softer, n = 1 polytropic 
EOS. 



